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ABSTRACT 



Context. Stochastic acceleration is thought to be a key mechanism in the energization of solar flare electrons. 

Aims. We study whether stochastic acceleration can reproduce the observed soft-hard-soft evolution of the spectral features of the hard X-ray 
emitted by suprathermal electron. We pay special attention to the effects of particle trapping and escape. 

Methods. The Fokker-Planck equation for the electron distribution is integrated numerically using the coefficients derived by Miller et al. 
for transit-time damping acceleration. The electron spectra are then converted to photon spectra for comparison with RHESSI observation of 
looptop sources. 

Results. The presence of particle escape softens the model spectra computed in the stochastic acceleration framework. The ratio between the 
efficiency of trapping and acceleration controls the spectral evolution which follows a soft-hard-soft pattern. Furthermore, a pivot point (that is, 
a common crossing point of the accelerated particle spectra at different times) is found at around 10 keV. It can be brought into agreement with 
the observed value of 20 keV by enhanced trapping through an electric potential. 

Conclusions. The model proposed here accounts for the key features observed in the spectral evolution of hard X-ray emission from looptop 
sources. 

Key words. Sun: flares ~ Sun: X-rays, gamma rays - Acceleration of particles 



1. Introduction 

Solar flares emit large amounts of continuum hard X-ray 
bremsstrahlung from energetic electrons up to several hundreds 
of keV. The thermal energy of the ambient electrons in the 
corona is substantially lower, with electron temperatures ob- 
^ ^ served up to 6 MK in nonflaring active regions (Brosius et al. 
• , 119961 Benz & Grigis EUU^ . Clearly, a mechanism is needed to 
' accelerate the electrons out of the thermal distribution up to rel- 
ativistic energies. More than 10-" erg (Lin et al. 2003 Emslie 
- - ' et al. l2004l Kane et al. l20()5l of energy can be transferred into 
non-thermal electrons during a major solar flare over a period 
of a few hundred seconds. 

These numbers set stringent requirements to acceleration 
models, and it is by no means a simple task to identify the 
physical processes involved. Most of the proposed models fall 
into three broad classes: electric DC field acceleration, stochas- 
tic acceleration and shock acceleration (see e.g. the review by 
Aschwanden 2002 j . 
Miller et al. 



il996l hereafter MLM) proposed a stochastic 
acceleration mechanism where electrons are energized by small 
amplitude turbulent fast-mode waves, the transit-time damp- 
ing model. MLM showed that their model could successfully 
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account for the observed number and energy of electrons ac- 
celerated above 20 keV in subsecond spikes or energy release 
fragments in impulsive solar flares. However, they made no at- 
tempt to explain the observed hard X-ray spectra (which are 
softer than predicted by the transit-time damping model) and 
did not consider spectral evolution. The MLM approach does 
not account for escape and transport processes. 

What is the effect of escape on the electron spectrum in 
the accelerator? In stochastic acceleration, each electron de- 
scribes a random walk in energy space. The effect of escape 
is to shorten the dwell time of the particles in the accelerator, 
and therefore it limits the average energy of the electrons. This 
results in a softer spectrum. Therefore, we need to take in ac- 
count escape to compare quantitatively the model predictions 
with hard X-ray observations. Furthermore, transport effects 
account for the modification of the electron spectrum from the 
time they leave the accelerator until they reach the hard X-ray 
emitting regions. The usual interpretation is that the particles 
are accelerated near the top of magnetic loops, but most of the 
X-rays are emitted in the loop footpoints, where the density is 
much larger. Nevertheless, looptop sources are also observed 
(Masuda et al. '1994l. 

Observations from the Reuven Ramaty High Energy 
Imaging Spectrometer (RHESSI, Lin et al. 2002 1 deliver hard 
X-ray spectra of both footpoint and looptop sources (Emslie 
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Table 1. Observational constraints for looptop sources (taken 
from Battaglia & Benz 12006 1 



Parameter Description Average Range 



Pivot-point energy e. 


20 


keV 


16- 


-24 


keV 


Pivot-point flux" J, 


2 




1-4 




Mean temperature 


22 




18- 


-25 


MK 


Nonthermal fitting range: 












Lower energy 


25 


keV 


20- 


-30 


keV 


Upper energy e^ax 


60 


keV 


40- 


-80 


keV 



In units of photons cm ^ s ' keV 



et al. l2003l Liu et al. l2004t . Looptop sources are better suited 
for comparison with accelerator models than footpoint sources, 
since one does not need to take into account the spatial trans- 
port from the accelerator The emission from the modeled elec- 
tron distribution can be directly compared with the observed 
spectra of the looptop source. 

BattagUa & Benz J2006t present a careful and comprehen- 
sive study of the spectral evolution of looptop and footpoint 
sources for 5 well observed RHESSI flares. They find that loop- 
top sources systematically show the soft-hard-soft (SHS) be- 
havior observed spectroscopically for the total flare emission 
(Grigis & Benz .2004. and references therein). On the other 
hand, some footpoint sources do not follow the SHS pattern. 
These observations suggest that the SHS behaviour is a prop- 
erty of the acceleration process. 

The observed SHS behavior manifests itself as a tight cor- 
relation between the curves describing the time evolution of the 
negative spectral index y and the normalization of the non- 
thermal power-law component of the photon spectra, measured 
at a fixed reference energy Eq. RHESSI observations show a 
linear relation between y and log Je^ (Grigis & Benz .2004 j . 
This is equivalent (Grigis & Benz I2005I I to the presence of a 
pivot point, that is, a common point where the power-law com- 
ponent of the photon spectra observed at different times inter- 
sect. Therefore, the position of the pivot point can be found 
observationally by fitting the regression line to y and log Je„ ■ 
Since the observations show that during the impulsive phase of 
most flares the scatter around such a model is relatively small, 
such fittings provide us with two new meaningful observational 
parameters which can be used for comparison with theories. 
The best-fit pivot-point energy can be measured by RHESSI 
with an accuracy of 10-20% (Battaglia & Benz"2006'). 

The measured energies of the pivot points lie in the range 
of 18-24 keV for looptop sources and 13-15 keV for foot- 
point sources. In Tabled we summarize the observational re- 
sults from Battaglia & Benz QQQ6J which we will use to con- 
strain the transit-time damping acceleration model. 

In this paper we study how well the transit-time damping 
model endowed with an escape mechanism can account for 
the hard X-ray spectra of the looptop sources seen by RHESSI 
and, in particular, for their spectral evolution. We investigate 
whether the spectrum of the photons emitted by the modeled 
electrons in the accelerator shows SHS behavior and how well 
this quantitatively agrees with the observations. 



We proceed as follows: a modified version of the MLM 
model with a better characterization of electron trapping is pre- 
sented in Sect. |3] The evolution of the electron energy distri- 
bution function is given by a diffusion equation which is in- 
tegrated numerically and transformed into photon spectra. In 
Sect. 13 we present the resulting photon spectra and compare 
their behaviour with the observations. The results are discussed 
in Sect. Island conclusions are drawn in Sect.|6l 

2. Pivot-point theory 

In this section the mathemathical foundations needed for the 
comparison of the model results with the observations are pre- 
sented. In particular we explain what are the consequences of 
the correlation in time between flux and spectral index. By the 
introduction of a new parameter, the pivot-point, the time de- 
pendence can be eliminated, thus simplifying the analysis. 

The theory is presented in a general context, and can be 
applied to any kind of spectra (for example, photons or elec- 
trons) described by a power-law in some energy range. The 
starting point of this analysis is the presence of a correlation 
between the logarithm of the flux, log , measured at energy 
E() and the spectral index 6, where the time dependent spectrum 
is given by 

F(E,t) ^FE„(t) i—j , (1) 

for f in a time interval Tint (for instance, the duration of a flare 
emission peak). Panel (a) of Fig.^shows the ideal case where 
the correlation in time between 6 (lower curve) and log Fe„ (up- 
per curve) is perfect: the two curves are parallel. In such a case, 
a plot of 6 vs. log Feo is linear. This shown in panel (c). The 
relation between these two parameters can be expressed as 

6(t)^a-\ogFE,(t) + h, (2) 

where a 0, and the parameters a and b do not depend on time. 

Some points with selected values of 6 are marked by spe- 
cial symbols in Fig. [2 circles for 6-2, triangles for 6 = 4, 
squares for 6 = 6 and stars for <5 = 8. In panel (b) the spectra 
corresponding to these values of the spectral index are plot- 
ted together It is evident that all these spectra intersect in a 
common point, the pivot point. The coordinates of the pivot 
point are defined as (£», F^). The presence of a linear relation 
between 6(t) and log FE„{t) is equivalent to the presence of a 
pivot-point. This means that panel (a) in Fig.^implies both (b) 
and (c). A detailed proof follows. 

2.1. Proof of the equivalence between pivot point and 
correlation of 6 with log F 

The following statements are equivalent: 

(i) There is an Eq such that F{Et, fi) = F(Et, t2) =: F^ 
for all fi, f2 e Tint- 

(ii) There exist constants a i^Q and b such that 

6(t) = a\ogFE,(t) + b (3) 
for all t e Tint- 
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Fig. 1. Schematic illustration of the basic properties of a spectral evolution. Panel (a): time evolution of the logarithm of the 
flux measured at reference energy £0 (upper curve) and of the spectral index (lower curve) in the case where they are perfectly 
correlated. Panel (b): Spectra for different values of the spectral index, crossing at the pivot point. Panel (c): Linar dependence 
of the spectral index on the logarithm of the flux. In all panels, the stars, squares, triangles and circles mark points with spectral 
indices of, respectively, 8, 6, 4 and 2. 



The parameters a and b are given by 
1 



log(£'. /£■()) 

- log 
logiEJEo) 



We give first the proof that (i) 
for all f € Tint, Eq. Q yields 



logF. = log 



= log Fe„ - SlogiEJEo) 



which can be solved for 6, yielding 
log Feo logF. 



6 = 



log(£',/£'o) log(£',/£'o) 
= a log FE,+b, 



where a for finite values of E»,Eq. 

Now we prove (ii) => (i). Let us define - E()exp(l/a) + 
(•4-) £■(). Then using Eqs. @ and Q we get 



(5) log ^(^•) = log Fe,-5- log y 1^ 



(ii). Because F{Et,t) - 



(7) 



log Fe„ - - 
a 

b 

a ' 



(9) 

(10) 
(11) 



(6) which does not depend on 6 or Fe^- Therefore, (£,, F,) are the 
coordinates of the pivot point. Q.E.D. 



2.2. Implications 



We have established above formally the equivalence between 
(8) the presence of a pivot point and the linear relation between the 
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Energy 

Fig. 2. Spectra (mock data, continuous lines) turning over at 
low energy. The power-law extensions of the spectra (dashed 
lines) intersect in a virtual pivot point, with larger flux than 
the one from the spectra. The dashed-dotted line represent an 
additional thermal component, effectively hiding the spectral 
behavior at low energies. 



We emphasize that the presence of a pivot-point is a local 
property of a family of curves at energy £0, in the sense that 
it can be seen as a convergence point for the tangents in Eq 
to the family of curves (in logarithmic space). For power-law 
functions, the coordinates of the pivot point are independent of 
the energy Eq, because their logarithmic derivative is constant 
as a function of energy. Other functions exist for which there is 
a different pivot point for each value of Eq. 

For the sake of completeness, we report here the differential 
equation whose general solution is the family of curves yield- 
ing a pivot point (in the sense of a convergence point for the 
tangents in linear space) in x*(x),y*{x): 

^+ ^ ^ ^'^"^ , (12) 
dx x*(x) — X x*(x) — X 

where in our case, x - \ogE, x* - \ogE*, y - \ogF and 
/ = logF*. 

In the special case x* - const, and y* - const, (that is, a 
unique pivot point for all energies), the general solution of Eq. 
(I12> is given by y - y* - m - (x- x*), the family of all straight 
lines intersecting at x* and y* . 

3. The Model 

The transit-time damping model is described in great detail by 
MLM. Here we first summarize the important features of the 
model presented by MLM, and then describe our modification 
of the model. 



spectral index and the logarithm of the flux measured at some 
fixed energy. We note that the shape of the curves describing 
the time evolution of 5 and Feq does not influence the results 
described above. The time dependence can be eliminated: this 
simplifies the analysis and allows the new parameter, the pivot- 
point position, to be used in the comparison of observations 
with theoretical models. Conversely, if the observation show 
that 6 and logF^u do not correlate well in time, we will not 
find any pivot-point like behavior in the spectral evolution. 

Our derivation assumes that the spectra are represented 
by power-laws at all energies. In this case, at energies lower 
than the pivot-point energy, the linear relation between flux 
and spectral index changes sign and higher fluxes would cor- 
respond to steeper spectra instead of harder ones. However, the 
observed spectra are not given by power-laws at low energies, 
because steep spectra are not integrable and therefore would 
imply an infinite number of particles, which is not physically 
possible. Therefore a deviation from the power-law form is ex- 
pected in the form of a turnover at low energies (Saint-Hilaire 
& Benz l2003i Sui et al. l2005t . This means that it is well pos- 
sible that the spectra do not extend all the way to the pivot 
point energy, even if a correlation between the spectral index 
and the flux is observed at higher energies. In such a situation 
the pivot point is virtual, in the sense that its flux is higher than 
the spectral flux at the same energy. In this case the analysis 
given above retain its validity, but the spectra used should be 
understood as the power-law extensions of the physical spec- 
tra. Such a situation is represented in Fig.|2l 



3.1. The Transit-Time Damping Model 

The source of the flare energy is an event of magnetic re- 
connection, where magnetic field energy is transferred, by an 
unspecified mechanism, into small-amplitude magnetoacoustic 
fast mode waves. These waves have very large wavelengths, 
comparable with the length scale of the global magnetic field 
restructuring induced by the reconnection process. 

A cascading process then takes place, producing waves 
with larger wave vector k and isotropic distribution. This pro- 
cess transfers wave energy density from the low k into the large 
k waves. The electrons interact with the parallel component of 
the magnetic field of the large k waves. An electron of velocity 
V can efficiently exchange energy with the waves if the reso- 
nance condition A:||V|| - ll) - kvp, is satisfied, where the subscript 
II denotes the component parallel to the ambient magnetic field, 
(jj is the wave frequency and va is the Alfven speed. The same 
condition can also be expressed as yu/yv = va, where yu and 77 
are the cosines of the angle between the ambient magnetic field 
and the direction of propagation of the electrons and the waves, 
respectively. From this form of the equation it is clear that only 
electrons faster then v'a can be accelerated by this mechanism. 
The perpendicular velocity Vj^ is not affected by these interac- 
tions. 

Since the acceleration only occurs in the direction of 
the ambient magnetic field, the electron distribution becomes 
strongly elongated in parallel direction in the absence of pitch- 
angle scattering. This reduces considerably the efficiency of 
transit-time damping acceleration, because for fast particle 
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with V » v'A and \f4 « 1 the resonance condition can only be 
satisfied by waves with \ri\ < va/v, which represent only a small 
fraction of all the available phase space. The acceleration effi- 
ciency can therefore be increased by the presence of a mecha- 
nism which isotropizes the electron distribution, allowing parti- 
cles with small |yL(| to interact with waves with \ri\ > va/v. MLM 
assume that an isotropization mechanism is present, but do not 
address the details of the process. 

The model described by MLM is self consistent in the sense 
that it describes the evolution in time of both the electrons 
(which are accelerated) and the waves (which are damped). The 
efficiency of the electron acceleration by the waves is deter- 
mined by the total energy density of the waves Ut and by the 
shape of the wave spectrum: large k waves exchange energy 
with the electrons faster because the transit time between the 
magnetic perturbations gets shorter. The coefficients describ- 
ing the diffusion of the electrons in energy space depend on the 
product of Uj and {k}, where the latter represents the average 
wave-vector of the waves. 

The evolution of the electrons is described by the Fokker- 
Planck equation for the electron distribution in energy space. 
MLM consider the effects of the acceleration by the waves and 
of collisions with the ambient plasma. 

3.2. Transit-Time Damping and Escape 

In this paper, we study the electron and photon spectra pro- 
duced by a simplified transit-time damping acceleration model. 
We do not address here the problem of the time it takes for the 
accelerated electron distribution to reach equilibrium, consider- 
ing that MLM have established this time to be sufficiently short. 
Therefore we need not to account explicitely for the MHD cas- 
cading process, but we will instead keep Uj-ik) as a free param- 
eter of our model. This means that only one partial differential 
equation describes the time-evolution of the system instead of 
a set of two coupled equations. 

A short description of the physical parameters upon which 
the model depends is given in Table |2l where a default refer- 
ence value for each parameter is listed. In the following these 
default values will be implied if different values are not ex- 
plicitely mentioned. 

The transit-time damping acceleration model describes the 
evolution of the electron population in energy and time. We de- 
note by N{E) dE the number of electrons per cubic centimeter 
with energy in the interval dE around E, where E is the di- 
mensionless particle energy in units of m^-c^. The total electron 

/>oo 

density is n = J^^ N(E) dE. The evolution in time of the dis- 
tribution function N(E, t) is given by the convective-diffusive 
equation 



Table 2. Description of the parameters used in the model. 



dN 



1 



(DCOLL + Dt)N 



dt 2 5£2 

-S(E)-N + Q(E) 



d_ 

dE 



(AcOLL +At)A^ 



(13) 
(14) 



The coefficients Dj and At describe the diffusion and con- 
vection in energy space due to the interaction of the electrons 



Param. 


Description 


Default value 


B 


Ambient magnetic field strength 


500 G 


n 


electron density 


10'" cm ^ 


T 


Ambient plasma temperature 


10 MK 


log A 


Coulomb logarithm 


18 




Proton gyrofrequency 


4.78 MHz 


Tn 


Time unit Q.^^ 


2.09 • 10"' s 


va 


Alfven speed 


0.034c 


Uj 


Energy density of the accelerating 






waves 


2 erg cm"^ 


Ub 


Energy density of the ambient 






field (|fi-/^) 


10"* erg cm"' 


(k) 


Average wavenumber of the 






turbulence 


4.8 • 10-" cm-' 


Iacc 

T 


Acceleration parameter: 
Ut cik) 

Escape time 


(free) 
(free) 



with the waves. The coefficients Dcoll and Acoll represent 
the effect that collisions against the ambient plasma have on 
the distribution function. We include the effects of the escape 
by a sink term S (E) and a source term Q{E) described below. 

The acceleration coefficients are given by MLM, and we 
report them here written in dimensionless form: 



At(£) 



-/3llACcy/3g(J3) 



Dj{E) = -y6i/ACcr¥/08) 



where 

r = 1 +£■, 



y8= a/^ 



r 



m = -1.25 -(1+2^2) log ^ + ^2 ^0.25^ 
1 



g(J3) = —(Aflog^-t + l) + m 



(15) 
(16) 

(17) 

(18) 
(19) 



The acceleration coefficients vanish for E < E/^, where ^a is 
the kinetic energy of a particle with speed /Sa - va/c, because 
the resonance condition cannot be satisfied below the Alfven 
speed. 

The coefficients for Coulomb collisions of the acceler- 
ated particles with the thermal background population (Spitzer 
1962 Trubnikov 1965 1 are given by : 



Acoll = - v n [i//{x) - i/f' (x)] E'"^ 

r, o ^1/2 
£>C0LL = 2 vn E ' 

X 

with 

ne log A Tu 



kBT ' 



Jo 



dt . 



(20) 
(21) 

(22) 
(23) 



The energy dependence of the acceleration and coUisional co- 
efficients is shown in Fig.|3] 
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Fig. 3. Upper panel: the absolute value of the dimension- 
less Coulomb (subscript coll) and transit-time damping (sub- 
script t) convection coefficients. Lower panel: the dimension- 
less Coulomb and transit-time damping diffusion coefficients. 
The arrows mark the thermal energy of the electrons in the 
plasma and the kinetic energy of electrons moving with the 
Alfven speed v^. 

The sink term is given by: 
S(E) = (24) 

T 

where t is the escape time. This represents a purely kine- 
matic escape for a spatially homogeneous electron distribu- 
tion (leaky-box escape). We can understand this term in the 
following way: let's assume that the plasma is contained in a 
cylinder with length L and cross section A, having a spatially 
homogeneous electron density n with a total number of parti- 
cles A^roT = J n dV - n ■ L ■ A all having the same speed v. 
After a time interval dt, the number of particles which have 
escaped the cylinder (assuming free streaming without colli- 
sions) is dNjcfT - n-v- A - dt. The rate of change of the density 
K dnjdt - n ■ vjL. Defining the escape time t - Ljc yields 
Eq. J24t . where the factor Th = O.^ comes from the transfor- 
mation of the equation into dimensionless form. 

In a more general case, where collisions have an effect, the 
particles cannot stream freely out of the accelerator, and they 
will need a longer time to escape. In this case, the escape time t 
can be larger than Ljc and the energy dependence on the colli- 
sion frequencies will cause t to vary with the energy. Since we 
do not address here the physical details of the escape we will 
assume for now that t{E) is constant for all energies. 
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Fig. 4. Accelerated electron distributions with different values 
of the power-law index resulting from changes in Ij - Iacc ■ t. 
The dashed curve represents the ambient Maxwellian distribu- 
tion. The two dotted lines indicate the energy range used for 
the computation of the power-law index 6 shown above each 
spectrum. 



Under the influence of collisions, we are not allowed to use 
L as the length of the path traversed by a particle before it can 
leave the accelerator: L(E) = ct(E) denotes the effective length 
for particle escape. It is best to think of the escape time r as a 
general parameter describing the strength of the particle trap- 
ping in the accelerator: the trapping becomes the more efficient, 
the larger t. 

The source term QiE) is needed to keep the number of par- 
ticles in the accelerator constant. Physically, these are electrons 
supplied by a return current, needed to keep electrically neu- 
tral the region where the electrons are accelerated. We assume 
that the electrons supplied by the return current mechanism are 
"cool" in the sense that their distribution is comparable with 
a thermal spectrum and with the temperature of the ambient 
plasma. Therefore we have 



e = no ■ Nmb(E) , 
where 



(25) 



(26) 



is the Maxwell-Boltzmann distribution normalized to unity, 
and ho = J S N dE is tiie rate of escaping particles. 

3.3. Method of Solution 

Equation J13I I is solved numerically using the Crank-Nicolson 
finite differences scheme in logarithmic energy space, which is 
well suited for diffusion problems. This method is the same as 
used in MLM and gives accurate results with fixed steps in time 
and log E. 
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3.4. From Electrons to Photons 

The numerical solution described above yields electron spectra. 
For a meaningful comparison with the observations we need 
photon spectra. They are generated by bremsstrahlung from the 
electron distribution. Since the model spectra are in equilibrium 
(particle losses by escape and coulomb collisions are compen- 
sated by the acceleration) a thin-target emission is computed. 

We convert the energy differential electron flux distribution 
F{E) in electrons cm"^ s ' keV ' into the photon spectrum 
J{s) observed at Earth in photons cm"^ s ' keV ' using the 
equation for the thin-target bremsstrahlung emission: 

7(e) de=-^ N{E)do-{E,E)dE , (27) 

where V is the source volume and R is the distance from the 
Sun. The cross section dcr(s,E) used is the fully relativis- 
tic, spatially integrated Bethe-Heitler formula (Bethe & Heitler 
I1934> without further approximations. This should adequately 
represent the emission process of the computed electron spec- 
trum in a looptop source. 

We will consistently use the notation 6, j, F^^, J^.^, 
and 7» for, respectively, the electron spectral index, the pho- 
ton spectral index, the electron density at energy Eq, the photon 
flux at energy sq, the electron pivot-point energy, the photon 
pivot-point energy, the electron density at and the photon 
flux at e». 

4. Results 

Section 14.11 presents the detailed properties of the numerical 
solution of Eq. (I13> for the standard set of parameters ("default 
values" given in Table |2j and the energy independent escape 
model. In a second step (Sect. l4~2t we will proceed to study the 
dependence of the results on the values of the model parameter 
and, finally, examine an alternative escape model (Sect. l43t . 

4.^. Results for the Default Values of the Model 
Parameters 

In this section the default values for the model parameters pre- 
sented in Tableware used. In this case, t and /acc are the only 
free parameters of the model, and the acceleration coefficients 
are much larger than the Coulomb coUisional coefficients in the 
energy range above 10 keV (see Fig.|3. Therefore, we expect 
that the electron spectrum in this energy range depends only 
on the value of the product Ij - I acc ■ t (cf. Eq. I13> . and we 
can basically work with one free parameter. For some choices 
of the model parameters, this approximation will not be valid 
anymore, so we will relax this assumption in Sect. 14.21 

Figurel^presents the electron spectra resulting from the nu- 
merical solution of Eq. i\3l for 4 different values of which 
yield different values of the electron spectral index 6 fitted in 
the energy range 30-80 keV. The harder spectra are the ones re- 
sulting from a higher value of It- The dependence of the spec- 
tral index on It is shown in the top panel of Fig.|5l Note that 
in the case where acceleration is weak or escape strong (that is. 



7 

It is small) very soft spectra are produced. For increasing trap- 
ping and acceleration efficiency (It larger) the spectra become 
harder The middle panel of Fig. |5] shows the dependency of 
the electron flux at 50 keV, F^q, on 7^. Here the flux is larger 
for larger It- During a flare. It will change as more energy is in- 
jected into turbulence waves, and therefore both 6 and F^o will 
change. This explains qualitatively the soft-hard-soft effect: as 
It increases, 6 will decrease and F^q will increase. It is shown 
in the bottom panel of Fig.|5] where 6 vs. Fso are plotted as a 
function of the parameter 7^. 

Since It contains the physics of acceleration and escape, 
the time evolution of the spectral index and flux in our model 
is a direct consequence of the changes in 7^. For example, if the 
energy density of the turbulent waves grows, reach a maximum 
value and then decreases, the electron spectrum will harden un- 
til peak time and soften again. 

To check the quality of the numerical solution, we obtain 
an approximated analytical expression for the function 6(It), 
whose computation is explained in Appendix IaI The approxi- 
mate curve is plotted in the top panel of Fig.|5las a dashed line 
and agrees well with the numerical solution. Because of the lin- 
ear nature of Eq. il3\ . the method exploited in App.lAl cannot 
be used to get an approximate solution for Fsq(It). 

The dashed line in the bottom panel of Fig.|5]represents the 
best fit of a logarithmic function (linear in log-lin space) to the 
data for 6 in the range 2-8. Such a relation between the electron 
flux at a reference energy and the spectral index indicates the 
presence of a pivot point in the electron spectra, as shown in 
section|2 From the slope and normalization of the dashed line 
we can get values of the pivot-point energy E, and flux 7"* for 
the electron distributions. 

More precisely, if 

6^alogFE,+b, (28) 
then 

= Eoe^'" and F, = e"*/° (29) 

as follows from Eq. 0. The dashed line in Fig. |5j bottom, 
corresponds to E, - 7.65 keV and 7^* = 2-10** elec- 
trons cm ke V " ' . 

For a meaningful comparison of these results with the ob- 
servations, a value for the pivot point of the photon spectra is 
needed. 7, and e, can be computed by fitting the regression line 
of y vs. log 7e„, in the same way as and 7^*. Note that while 
Jt depends on the source volume, e» can be immediately com- 
pared with the observations. Using the photon spectral index 
computed in the energy range 30-50 keV, the regression (re- 
stricting y to the range 3-9) yields the pivot-point coordinates 
£, = 5.0 keV and 7, = 1.9 ■ 10^^ ■ Vn photons cm^^ s"' keV"', 
where ¥21 is the volume of the source in units of 10^^ cm^. 

The above value of the pivot-point energy is too low com- 
pared to the looptop source observations by Battaglia & Benz 
I2OO6 1 summarized in Table [2 reporting e» x 20 keV. This 
means that if the photon spectral index of our model varies 
from, say, 3 to 6, the non-thermal flux is subjected to an ex- 
cursion which is too large to account for the observed behavior 
of the spectrum in solar flares. 
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Two possible solutions of this problem exist: we can ex- 
plore the parameter space to find out if there is some combi- 
nation of the parameter values which indeed produces a higher 
pivot-point energy e,, or we can try to modify the escape term. 

4.2. Exploration of the Parameter Space 

The set of default parameters used in the previous section yields 
a pivot-point energy e, which is lower than the observed value. 
However, the pivot-point energy depends on model parame- 
ters like T, n etc. Therefore a better coverage of parameter 
space than the example reported above is needed to assess the 
range of variability of £». The pivot-point energy depends on 
the temperature of the ambient plasma T, its electron density 
n and its magnetic field strength Bq (indirectly, through the 
Alfven speed). We therefore compute the equilibrium solutions 
of Eq. J13> for a large number of combinations of different val- 
ues of these parameters, each set yielding a value for e». In this 
way we sample the function e,(r, «, ■ ■ ■)• 

There is however a complication: for some choices of the 
parameters, the approximation that the spectral index and flux 
depend only on = r ■ /acc will not be valid anymore. This 
approximation becomes invalid if AcollMt or Z)coll/£'t ap- 
proach unity. We note that both ratios are roughly proportional 
to rP' at high energies, and therefore the approximation be- 
comes bad for large values of the density. In this case, the es- 
cape time T and the acceleration parameter /acc separately in- 
fluence the spectral index and flux, and the pivot-point energy 
will depend on the t vs. /acc relation. Since we have not mod- 
eled the physical mechanism responsible for the trapping of 
particles, the details of how r depends on /acc are not known. 

We estimate the range of pivot points by choosing several 
representative paths in the r, /Acc-plane (as shown in Fig. |6|l 
and computing the pivot-point coordinates resulting from spec- 
tra computed along these paths. 

Recalling that in first approximation the spectral index de- 
pends only on the product Ij - /acc ■ t, the paths are required 
to start from a point on the line r ■ /acc = 10"^ and end in a 
point on the line t ■ /acc - 10 The region between these 
two delimitation lines is shown in gray in Fig.|6l These limits 
are chosen to ensure that the computed photon spectral indices 
fall into the observed range. For the standard values of the pa- 
rameters, the region above the gray area yields spectral indices 
harder than about 1.5 and the region below yields y softer than 
about 15. 

There are however more constraints on the values of /acc- 
This parameter is the product of the dimensionless wave num- 
ber and the energy density in the turbulent fast mode waves. In 
our simplified version of transit-time damping acceleration, the 
time evolution of the spectral energy distribution of the waves 
is not computed. Therefore we set the maximum value of /acc 
= 6 ■ 10""*. It is the maximum value that is obtained 



at / 



ACC 

by MLM during the cascading process (read from the plot in 
Fig. 6b in MLM). Since /acc controls the acceleration effi- 
ciency, values which are much smaller then yield a weak 
and slow acceleration, contrary to the observations. Larger val- 
ues of /acc could result from a higher energy density in the 



waves, but in that regime the wave amplitude becomes so large 
that nonlinear effects are likely to a lead to a breakdown of the 
physical model used by MLM to compute the acceleration co- 
efficients. Thus values of /acc much larger than I^qq do not 
necessarily improve the acceleration efficiency. 

For these reasons we additionally require that all the paths 
end on the point {t^^^, ^acc^)' where t"^^^ is the point on the 
upper delimitation line with ordinate I^qq- Furthermore, we 
also require the paths to be monotonically increasing functions 
of both coordinates (this guarantees that the product Ij is also 
monotonically increasing). Our reference paths are thus given 
by 



/acc - /acc 



\^MAX j 



(30) 



for different values of a between and oo. They represent 
straight lines in the logarithmic /acc vs. t plot shown in Fig.|6l 



intersecting in the point with coordinates (r 



MAX rMAX 
^ACC 



The results are shown in Fig. where the value of the 
pivot-point energy e» is plotted for 3 diff'erent values of the 
plasma temperature (T - 10,20,30 MK) and several differ- 
ent values of the plasma density between n = 10^ cm and n - 
10^' cm The different points for each temperature and den- 
sity represent the range of variation of e» for different values of 
a as explained above. We have used a -Q, V2- 1,1, V2+ 1 , oo. 

The pivot-point energy increases with temperature. The 
density dependence of the pivot-point energy is weak, while 
the pivot-point flux decreases at higher densities. The num- 
ber of accelerated particles becomes very small at densities 
of 10'^ cm"^ and higher, because the particles lose energy in 
the collisions faster than they gain energy from the waves. The 
spread of the pivot-point energies shown in Fig. increases 
at higher densities, as expected from the high-density break- 
down of the approximation that the pivot point only depends 
on /acc ■ T. 

4.3. Alternative Escape Modeling 

As an alternative way for explaining the high values of e,, we 
now modify the escape mechanism. For instance, if an electric 
potential Ve is present between the accelerator and the foot- 
points, electrons with kinetic energy lower than the threshold 
energy Ej = cVe will not be able to cross that barrier and will 
therefore not leave the accelerator. Such an electric field is ex- 
pected to exist and drive the return current of electrons from 
the chromosphere. In this scenario the trapping is ideal below 
Ej, and r becomes infinitely large. 

This means that the escape time is given by: 



T(£) = 



if E <Ej 
if E> Et 



(31) 



The results are shown in Fig.|8]for two values of Ej, 25 keV 
and 45 keV. Note that the value of the pivot-point energy e, in- 
creases with increasing Ej. The pivot point energies reach the 
observed values for Ej ^ 30 keV, although even at 45 keV 
some paths in the r, /Acc-plane deliver lower pivot-point ener- 
gies than observed. 
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Fig. 8. Same as in Fig.0 but with a modified escape term featuring perfect trapping below 25 keV (left panel) and below 40 keV 
(right panel). 



The points in Figs. and |S] with larger pivot-point energy 
correspond to the paths with the lowest values of a, where the 
escape time r changes faster than the acceleration parameter 
Iacc- 

A comparison of observed spectra of a looptop source (for 
the M3 flare of December 4, 2002) and the spectra computed by 
the stochastic acceleration model with Ej = 40 keV is shown 
in Fig. |9] The model parameters used were temperature T = 
20 MK and density n - 5 ■ 10"^ cm"-* and a - y/l - 1, yielding 
a pivot-point energy e» = 18.7 keV, similar to observed value 
of 18.1 keV reported by BattagUa & Benz (2006). The model 
spectra are represented by the dashed lines, the power-law fit 
to the data by the dotted line, and the sum of the model spectra 
with an isothermal component by the continuous line. 

In our model, the value of the total photon flux observed 
at earth depends linearly on the accelerator volume, which acts 
as a normalization factor for the model spectra. For the com- 
parison with the observed data, the volume can be freely cho- 
sen, but is assumed not to change along the path in the r, /acc- 
plane. In practice, that means that the volume can be chosen 
in order to match the data (or the power-law fittings) at, say, 
peak time, and automatically the spectra at all the times dur- 
ing the emission peak will match, provided that the pivot-point 
energies of the observed and the model spectra are the same. 

The low energy part of the observed spectrum shows a ther- 
mal emission much larger than the model spectra emission. 
This can be understood if particle acceleration takes place in a 
smaller volume than the one where the thermal emission takes 
place. For the events shown in Fig. |9] the filling factor of the 
accelerator amounts to around 10"^. 

The left panel of Fig.|9lshows a mismatch between the ob- 
served data and the model spectra below 30 keV. This is due 



to the effect of the escape model used: the suppression of elec- 
tron escape results in a hardening of the electron spectra below 
Ej - 40 keV. Particles with energy lower than Ej are able to 
extract more energy out of the waves until they are acceler- 
ated beyond Ej. A turnover occurs at around QJEj for photon 
spectra, due to the shape of the bremsstrahlung cross section 
do-(s, E). 

5. Discussion 

The important results are the following: 

- The transit-time damping mechanism with escape is able 
to provide photon spectra with spectral indices in the deka- 
keV energy range corresponding to the observed range for 
flares (y = 2 to 10, Dennis [T985I Battaglia et al. 20051 as 
the ratio between the trapping and accelerating efficiency 
is varied (Fig. 0}. This confirms that escape indeed does 
soften the spectrum for different values of densities and 
temperatures of the accelerator. 

- Assuming constant values of temperature and density in the 
accelerator (as may be the case during a hard X-ray peak 
with a duration of some tens of seconds) we are able to pro- 
duce a correlation between the photon spectral index and 
flux yielding an approximate pivot point in photon flux vs. 
energy space (Figs.|5]andE} 

- The energy of the pivot point in the simplest model is 
around 10 keV and is lower than the reported observed 
value of 20 keV. 

- An additional increase of the strength of trapping below 20- 
30 keV (as may be obtained from the presence of a potential 
barrier) increases the pivot-point energy to 20 keV (Fig.|8}. 
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Spectral Index: 6.3 Spectral Index: 3.1 




Energy (keV) Energy (keV) 

Fig. 9. Looptop source spectra observed during the M3 flare of December 4, 2002, represented by crosses (bin width and error 
bar, data courtesy of M. BattagHa). The dotted line is a power-law fit to the data (a low-energy turnover at around 15 keV) with 
spectral index y of 6.3 for the left panel and 3.1 for the right panel. In each panel, the dashed line is a stochastic acceleration model 
spectrum computed for model parameters: temperature T = 20 MK, density n - 5 - 10'*' cm"^, and threshold energy Ej - 40 keV 
(chosen such that the pivot-point energy s, matches the observed value of ~ 18 keV reported in Battaglia & Benz l20()6l . The 
continuous line is the sum of the model spectrum with an additional isothermal component, needed to fit the low-energy range of 
the spectra. 



- The potential barrier hardens the spectrum below its thresh- 
old energy. For large y, this flattens the spectrum more than 
observed in the range where the non-thermal emission still 
dominates the thermal emission (Fig.|5] left). 

- The observational upper limit for the energy of a turnover 
in the photon spectra lies around 15 keV (Saint-Hilaire & 
Benz 2005 1. Model spectra computed for threshold ener- 
gies below about 23 keV yield a turnover lower than 15 keV 
and pivot-point energies below 16 keV. 

What is the physical significance of the pivot point? We 
would like to point out first that there is not necessarily any: 
it can be thought of as a useful parameter which describes the 
variation in flux and spectral index of the observed spectra, and 
one has to explain this variation first in the context of a model 
such as done in this work. However, the fact that it appears at an 
energy which is comparable to the estimated turnover energy of 
the non-thermal component, as well as the energy below which 
the thermal component dominates is a suggestive fact. 

The model suggests that the main factor influencing the 
pivot-point energy is the number of electrons available at en- 
ergy comparable with it. We propose that, in general, the pres- 
ence of a large supply of electrons at some energy implies a 
similar energy for the pivot point, on the condition that the elec- 
tron collisional energy losses are negligible above that energy 
range, such that the acceleration process can eff'ectively ener- 
gize that population. It is not important if such a population 
exist at the beginning of the acceleration, but it must be present 



in the later stages and be maintained until the equilibrium is 
reached. 

6. Conclusions 

The model studied assumes that turbulent fast-mode waves are 
present in the accelerator region. The equilibrium spectrum 
of electrons subjected to Coulomb collisions and transit-time 
damping interactions with the waves are computed allowing 
for trapping/escape and replenishment of "cold" particles. As 
such the model does not distinguish between the accelerator 
and the radiator and considers the photon spectra emitted by 
the electrons in the accelerator itself. The relatively dense loop- 
top sources observed by RHESSI fit this scenario as acceler- 
ators/emitters, and the observed high densities support a sce- 
nario in which trapping plays an important role in acceleration. 
Therefore, spectral observations of these looptop sources may 
well deliver a snapshot of the electrons in the accelerator 

To match the low-energy part of the observed spectrum, we 
have to add an additional thermal component at the observed 
temperature. This component fills a much larger volume (factor 
of 1000). The heating may have occurred by earlier accelera- 
tion or may be the result of waves escaping from the accelerator 
region with energy density too low to yield significant energy 
gain per particle. 

The goals of this work were to test whether a stochastic ac- 
celeration mechanism can account for the observed soft-hard- 
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soft behavior and, in particular, to find a minimal set of mod- 
ifications to the transit-time damping mechanism yielding a 
pivot-point like behavior of the electron spectrum in the deka- 
keV range. Therefore we have not fully investigated the physics 
of escape, and in particular we leave open the question of the 
connection between the acceleration process and the trapping 
mechanism. 

The results shown in Section 0] confirm that the hardness 
of the spectrum of accelerated particles in a stochastic accel- 
eration model depends strongly on particle trapping. This can 
be understood by recalling that stochastic acceleration can be 
thought as a transport of particles in energy space due to a ran- 
dom walk process. While acceleration pushes particles toward 
higher energies, the escape acts against this flow, because fast 
particles are more likely to be lost, and replaced by particles at 
lower energies. Therefore, when the losses are stronger, the av- 
erage time a particle spends random-walking in energy space is 
shorter, thus gaining less energy. This results in softer spectra. 

We have used one specific kind of stochastic acceleration 
model and a simple escape scenario, but the physical mecha- 
nisms explored by the specific models are general enough that 
other stochastic models and escape terms are expected to fol- 
low the same general behavior leading to a soft-hard-soft effect, 
manifesting itself as a correlation of the time evolution of the 
spectral index and the non-thermal X-ray flux. 

A simple escape model cannot account for the observed 
value of the photon pivot-point energy, which is about a fac- 
tor 2 higher. The modification in the escape term featuring no 
escape below a threshold energy Ej ^ 40 keV can increase the 
energy of the pivot point to the observed values, but at the price 
of introducing a spectral hardening below 30 keV, which is not 
observed for large y. 

One possible interpretation of the threshold energy used in 
the alternative escape model is the presence of an electric field 
driving the return current. It was treated as a free parameter, but 
physically the escape process and the return current are linked 
by the electric conductivity of the loop legs. In this sense the 
model is not self consistent. The issues of transport of the es- 
caped particles to the footpoint and return currents are impor- 
tant (see for example Zharkova & Gordovsky '2005'^ and need 
to be addressed in future work to be able to derive footpoint 
spectra. 

Thus the model is not able to account quantitatively for 
all the observed features of the intricate spectral behavior of 
the hard X-ray emission from looptop sources of solar flares. 
Nevertheless, the model qualitatively accounts for the key fea- 
tures of the spectral evolution, going one step further toward the 
solution of the acceleration problem. This allows to interpret 
the impressive RHESSI observations showing spectral varia- 
tion in looptop sources. 

The simplifications introduced by the model may well be 
responsible for its shortcoming in reproducing the observed 
spectra for large y. Future improvements may include an 
isotropization process and magnetic trapping. This will require 
the extension of the model including at least one spatial dimen- 
sion and the pitch angle distribution of the particles. 
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Appendix A: Approximate analytical solution 

We compute here an approximate analytical solution in the en- 
ergy range kT/{mc^) « £ « 1 . In this range, we can neglect 
the influence of the Coulomb collisional coefficients Ac and Dq 
as well as the source term Q. In equilibrium Eq. il3\ becomes: 



(A.l) 



where the coefficients A, and D,- are the factors of order / occur- 
ring in the Taylor expansion of A and D around E. Their com- 
putation is tedious, but straightforward. Neglecting the term of 
order and setting y = 1 we get 



A) = K/3'{-5/4-log^) 
Di = 11 /4- 3 log ^) 

2K 



D2 



■ (-17/8 -3/2- log ^) 



(A.2) 
(A.3) 

(A.4) 

(A.5) 

(A.6) 



Ao = 2K/3i-l-\og^) 
2.K 

A, = y (-1/4- log ^) 

where |/32/Acc,^=y6A/ye. 

The solution N(E) will be approximatively given by a 
power-law around E - Eq 



N(E) = A^o 



(A.7) 



Plugging this function into Eq. iA.l\ yields a quadratic equa- 
tion for 6: 



C26 + ci6 + cq - 
with coefficients 



C2 



tDq 



Co = TD2 - tAi - JlE^ 



(A.8) 

(A.9) 

(A. 10) 
(A.ll) 



The solution of the quadratic equation yields 5 as a function 
of Ir = /acc ■ T and the model parameters. The asymptotic 
behavior of the solution is the following: 6 is constant for Ij — > 
00, and 6 is proportional to 1 / ^/Ti■ for — > 0. 
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Fig. 5. Top: The electron spectral index S (measured in the 
energy range 30-80 keV) and the photon spectral index y 
(measured in the energy range 30-50 keV) plotted against It- 
The dashed line represent the approximated analytical solution 
computed in Add. 1X1 

Middle: The electron density at 50 keV ^50 and photon flux at 
40 keV plotted against 7^. 

Bottom: 6 vs. F^q and y vs. 740 plotted parametrically as a func- 
tion of Ij. The dashed line are the best log-lin fit for 6 in the 
range 2-8 and for y in the range 3-9. 
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Fig. 6. Five paths (between filled circles) in the t, /acc- 
plane as described by Eq. J30l l. corresponding to arctana = 
0,;r/4,7!-/2, 37r/4,7r/2. The gray shaded region contains the 
points yielding photon spectra with a spectral index lying ap- 
proximatively in the 1.5-15 range. The two dashed lines have 
abscissa and ordinate of r'^^^ and /^cc^' respectively. 



Fig. 7. Dependance of the photon pivot-point energy e, on the 
electron density n for plasma temperatures T of 10 MK (stars), 
20 MK (circles) and 30 MK (triangles). The density values of 
the 10 and 30 MK data points have been shifted by 15% to 
avoid overlap with the other points. The spread arise from com- 
putation of the spectra along different paths as explained in the 
text. The gray shaded region represents the observed range of 
values of reported in Tabled 



